Onsager coefficients of a finite-time Carnot cycle 
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We study a finite-time Carnot cycle of a weakly interacting gas which we can regard as a nearly 
ideal gas in the limit of Th — T c — > where Th and T c are the temperatures of the hot and cold 
heat reservoirs, respectively. In this limit, we can assume that the cycle is working in the linear- 
response regime and can calculate the Onsager coefficients of this cycle analytically using the el- 
ementary molecular kinetic theory. We reveal that these Onsager coefficients satisfy the so-called 
tight-coupling condition and this fact explains why the efficiency at the maximal power ?? m ax of this 
cycle can attain the Curzon-Ahlborn efficiency from the viewpoint of the linear-response theory. 

PACS numbers: 05.70.Ln 
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I. INTRODUCTION 

Improving efficiency of heat engines has been a long 
challenge since the Industrial Revolution. Substantial 
progress of our knowledge came from Carnot 's insight: 
He conceived a mathematical model of an idealized heat 
engine, now called the Carnot cycle, and showed that 
there is an upper limit of the efficiency of all the existing 
heat engines which can be attained only when the engines 
are working infinitely slowly (quasistatic limit) to vanish 
the irreversibility. 

In the fundamental physics, properties of heat engines 
working at the maximal power have also been studied 
since the study by Curzon and Ahlborn [l|, Q (see also 
Q). Although the quasistatic Carnot cycle has the high- 
est efficiency, it outputs zero power because it takes infi- 
nite time to output a finite amount of work. By contrast, 
Curzon and Ahlborn considered a finite-time Carnot cy- 
cle which exchanges heat at a finite rate with the reser- 
voirs according to the linear time-independent Fourier 
law. Under the assumption of endoreversibility that irre- 
versible processes are occurred only through these heat 
exchanges, they derived a remarkable result: the effi- 
ciency at the maximal power r? max is given by 



(1) 



where Th and T c are the temperatures of the hot and cold 
heat reservoirs, respectively and the above 7/ca is usually 
called the Curzon-Ahlborn (CA) efficiency. 

Previously we studied a finite-time Carnot cycle of a 
weakly interacting gas which we can regard as a nearly 
ideal gas to confirm the validity of the CA efficiency from 
a more microscopic point of view [4, 5] . We performed ex- 
tensive molecular dynamics (MD) computer simulations 
of the finite-time Carnot cycle of the two-dimensional low 
dense hard-disc gas and measured the efficiency and the 
power for the first time. Our simulations revealed that 
our 7y m ax agrees with t/ca only in the limit of AT — > 
where AT = Th — T c , but exceeds r/cA at somewhat large 
AT. We also confirmed this behavior of r] max analytically 
using the elementary molecular kinetic theory. Therefore 
the phenomenological prediction of Eq. ([T]) seems to be 



valid only in the limit of the small temperature differ- 
ence for our finite-time Carnot cycle. Recently, Van den 
Broeck @ considered the heat engine described as the 
Onsager relations 



J\ — L\\X\ + L12X2, 
■h = L21X1 + L22X2, 



(2) 
(3) 



and showed that tjca is the upper limit of ?7 max in this 
heat engine (see also Sec. V). The recent studies @, H, 
B EH EI on the various theoretical heat engine models 
also support the results in [6]. These results seem to 
meet our previous result of the finite-time Carnot cycle, 
though it is unclear why our system realized the upper 
limit of r? m ax in AT — > 0. Moreover it is also unclear 
whether the finite-time Carnot cycle can be understood 
in the framework of the Onsager relations because the 
explicit calculations of the Onsager coefficients Lij for 
that cycle do not exist to our knowledge. 

In this paper, we apply the framework of the Onsager 
relations to our previous study of the finite-time Carnot 
cycle and analytically calculate the Onsager coefficients 
for it. Although there are a few analytic calculations of 
the Onsager coefficients for the steady state of heat en- 
gines such as Brownian motors 0, E2 , we believe that 
the present study is the first example of the calculation 
for the cyclic heat engine model where the two heat reser- 
voirs do not contact with the working substance simulta- 
neously. We will show that these Onsager coefficients sat- 
isfy the so-called tight-coupling condition and therefore 
we can give an explicit explanation why ?7 max of this cycle 
attains the CA efficiency in the limit of AT — > 0, as ob- 
served in 0, @ , from the viewpoint of the linear-response 
theory. We also perform the MD computer simulations 
to check the validity of our analytic calculations of the 
Onsager coefficients. 

The organization of this paper is as follows. First we 
introduce our finite-time Carnot cycle model in Sec. II 
and describe the molecular kinetic theory in Sec. III. The 
main result of this paper, the analytic calculations of 
the Onsager coefficients for our model, accompanied by 
the results of the MD simulations for checking the valid- 
ity of our analytic calculations, are shown in Sec. IV. In 
Sec. V, we introduce the general framework of the heat 
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FIG. 1: (a) Schematic illustration of a two-dimensional finite- 
time Carnot cycle model. Hard-disc particles as a weakly in- 
teracting nearly ideal gas are confined into the cylinder. The 
piston moves at a finite constant speed u and the thermal- 
izing wall with the length S is set on the left bottom of the 
cylinder during the isothermal processes, (b) Temperature- 
volume (T-V) diagram of a quasistatic Carnot cycle for a 
two-dimensional ideal gas. 



engine used in [6( and discuss the efficiency at the max- 
imal power of our finite-time Carnot cycle using those 
Onsager coefficients according to that framework. We 
summarize this study in Sec. VI. 



II. MODEL 

We first introduce a theoretical model for a finite-time 
Carnot cycle of a two-dimensional weakly interacting gas 
which we can regard as a nearly ideal gas 0, [f|. To 
mimic the weakly interacting nearly ideal gas, we confine 
a low dense N hard-disc particles with diameter d and 
mass m into a cylinder with rectangular geometry and 
let them collide with each other. (Fig. Q] (a)) The head 
of the cylinder is a piston and it moves back and forth 
at a constant speed u. The usual quasistatic Carnot cy- 
cle (u — > 0) consists of four processes (Fig. [U (b)): (A): 



isothermal expansion process (V\ — > V2) in contact with 
the hot reservoir at the temperature Th, (B): adiabatic 
expansion process (V2 — > V3), (C): isothermal compres- 
sion process (V3 —* V4) in contact with the cold reservoir 
at T c , (D): adiabatic compression process (V4 — * V\), 
where Vfe's (k = 1, • • • ,4) are the volumes of the cylin- 
der at which we switch each of the four processes. Vfc's 
satisfy the relations V 3 = (T h /T c )V 2 and V A = (T h /T c )Vi 
in the case of the two dimensional ideal gas. Since we 
regard our system of hard-disc particles as a nearly ideal 
gas, we apply these relations to our system. In the case 
of a finite-time cycle, we also switch each process at the 
same volume V% as in the quasistatic case. 

Defining (x,y) coordinates as in Fig. [T](a), we let the 
piston move along the ir-axis at a finite constant speed 
u. Here, we express the ir-length and the y-length of 
the cylinder as I and L, respectively and the volume of 
the cylinder as V — LI. Then, the re-length 1^ at the 
switching volume Vu can be defined as Ik = Vk/L. 

When a particle with the velocity v = (v x ,v y ) col- 
lides with the piston moving at the x-velocity ±u, its 
velocity changes to v' = (—v x ± 2u, v y ), assuming per- 
fectly elastic collision. Then the colliding particle gives 
the microscopic work m(|t>| 2 — |t>'| 2 )/2 = 2m(±uv x — u 2 ) 
against the piston. In the isothermal processes (A) and 
(C), we set the thermalizing wall with the length 5 at 
the position as in Fig. Q] (a). When a particle with the 
velocity v — (v x ,v v ) collides with the thermalizing wall, 
its velocity stochastically changes to the value governed 
by the distribution function 
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2n \k B Ti 



3/2 



v y exp 



2k B Ti 



(4) 



(—00 < v x < +oo,0 < Vy < +00, Ti (i = h in (A), c 
in (C))), where k B is the Boltzmann constant [13]. The 
microscopic heat flowing from the thermalizing wall can 
be calculated by the difference between the kinetic en- 
ergies before and after the collision. We sum up the 
above microscopic work and heat during one cycle. When 
Qh denotes the heat flown from the hot reservoir dur- 
ing (A) and Q c denotes the heat flown from the cold 
reservoir during (C), the total work W during one cy- 
cle is expressed as W = Qh + Q c - Then the power 
P and the efficiency rj can be defined as P = W and 
■q = W/Qh = P/Qh. Here, the dot denotes the value 
divided by one cycle period or the value per unit time 
throughout the paper. In our system, one cycle period is 
2(h ~h)/u. 



III. MOLECULAR KINETIC THEORY 

In this section, we review the results of the molecular 
kinetic theory of the finite-time Carnot cycle obtained in 
our previous work. The details of the derivation of the 
equations in this section are described in Q ■ 

If we assume that, even in a finite-time cycle, the gas 
relaxes to a uniform equilibrium state with a well-defined 
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temperature T sufficiently fast and the particle velocity 
v is governed by Maxwell-Boltzmann distribution at T, 
we can easily derive the time-evolution equation of T us- 
ing the elementary molecular kinetic theory. Such an 
assumption of the fast relaxation to a uniform equilib- 
rium state at a well-defined temperature T is valid when 
the energy equilibration in the cylinder due to the inter- 
particle collisions is much faster than the speed of the 
energy transfers through the thermalizing wall and the 
piston. This situation is surely realized when u is small 
and the interaction length S between the gas and the 
reservoir is sufficiently small. If we assume that the gas 
is sufficiently close to a two-dimensional ideal gas, the in- 
ternal energy of the gas can be approximated as Nk\ , / . 
Then, we can derive the time-evolution equation of the 
gas temperature T for each of the four processes (A)-(D) 



(A) : Nk B ^- = Oh ~ w c , (B) : Nk B ^- = -w e , 
at at 

(C) : Nk B ^- =q c - w c , (D) : Nka-jg = -w c , 



(5) 



where <?; = qi(t,T) (i = h in (A), c in (C)) is the heat 
flowing into the system per unit time in the isothermal 
processes and Wj — Wj(t,T) (j = e in the expansion 
processes (A) and (B), c in the compression processes (C) 
and (D)) is the work against the piston per unit time. 
By counting the number of the particles colliding with 
the thermalizing wall and the piston, we can derive the 
specific forms of qt and w c as 



3SNkBm-T) 2irk B T 
<li(t,T) = , t r / 1 \ \/ , (6) 



w e (t,T) = 



AirV(t) 



2muNL J A 2 T 



a T u 2 

A y- u +T 

V 7T Z 



dv x (^AVTv x 



(7) 



where A = yf2k B /m. w c is also obtained by changing 
u — > — u in Eq. (JT]). We can numerically solve Eqs. ([5]) 
for the entire cycle at various piston speeds. By using 
the final temperature of each process as the initial tem- 
perature of the next process repeatedly, we can obtain a 
steady cycle at a given u. In Fig. [2 we have shown the 
temperature-volume (T-V) diagram for the steady cycle 
at u = 10~ 4 (dotted line) and u = 10 -3 (dashed line), 
respectively. The solid line is the theoretical quasistatic 
line of a two-dimensional ideal gas. From this figure, 
we can see that as u becomes larger, the cycle deviates 
from the theoretical quasistatic line and the tempera- 
tures during the isothermal processes (A) and (C) relax 
to the steady temperatures Tjf (< T h ) and T c st (> T c ), re- 
spectively. If it is considered that the relaxation to the 
steady temperature is very fast, the isothermal process 
may approximately be divided into the relaxational part 




FIG. 2: Temperature-volume {T-V) diagram for the steady 
cycle obtained by numerically solving Eqs. ([5]) for u — 1CP 4 
(dotted line) and u = 10 -3 (dashed line). The solid line is 
a theoretical quasistatic Carnot cycle for a two-dimensional 
ideal gas. Since the dotted line is very close to the solid 
line, the case of u = 10 -4 may be regarded as the quasistatic 
cycle. The parameters used are N = 100, Th = 1, T c = 0.7, 
S = 0.05, m = 1, fe B = 1, h = 1, h = 1.5 and L = 1. 



and the steady part in the case of a finite-time cycle, 
where T(t) instantaneously changes from the initial tem- 
perature of the isothermal process to the steady temper- 
ature in the relaxational part and T(t) keeps the steady 
temperature in the steady part. These relaxational parts 
in the isothermal processes (A) and (C) are surely missed 
in the original model by Curzon and Ahlborn [l|, 0, EH • 
Moreover we note that the heat flow per unit time qi as 
shown in Eq. ([5]) is time-dependent even during T = 
through the volume of the cylinder after the relaxational 
processes in (A) and (C). Therefore the assumption of the 
time-independent heat flows in the original phenomeno- 
logical model in [l[ does not seem to be applied to our 
model. 

Although it may be impossible to find the exact solu- 
tions of Eqs. ([5]), we can find the analytic expressions of 
the steady temperature T? during the isothermal pro- 
cesses (A) and (C) as the solutions of dT/dt = by ex- 
panding as T?* = Ti + aju + a 2 u 2 + 0(u 3 ), assuming that 
u is small. These expansion coefficients aj , a 2 , etc. can 
be determined order by order. If we assume that the re- 
laxational process to is sufficiently fast, the heat 
flowing into the system during the steady part T(t) = TJ* 
can be calculated up to 0(u) as: 

Kh-h)/u 

Qh = / <lh(t,T?) dt 

Jo 

= Ql s -2mNA^(± + -^u\n^ (8) 

where the quasistatic heat for an ideal gas in the isother- 
mal process (A) is defined as = Nk B Th \n{y%/Vi). 
<3c* during T(t) = can also be obtained by chang- 
ing u -► -u, T h -> T c and V 2 /Vi -> V4/V3 in Eq. ©. If 
we neglect the contributions of the heat transfers through 
the thermalizing wall during the relaxational parts in (A) 
and (C), the steady part of the work during one cycle 
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W s 



Qh + Qf becomes 



and their conjugate fluxes as 



W st = Nk B (T h - T c ) In ^ - 2mNA ^ ( - 

Vl \TT 



L 
3S 



(9) 



If we assume that the adiabatic processes (B) and (D) 
even at a finite u satisfy the quasistatic adiabatic relation 
TV = const of a two-dimensional ideal gas, the additional 
heat transfers Q^f d and Q^ dd during the relaxational part 
in (A) and (C), respectively, can be estimated as 



-)add 



= -Nk 




(10) 



up to 0(u) and Q^ dd can also be obtained by changing 
u -c -u and T h <-> T c in Eq. ([TO]), Therefore the total 
Q h) Q c and IT become Q h = Q£+Q| dd , Q c = Qf+Qf d 
and W = W st + Q£ dd + Q* dd , respectively. 



IV. CALCULATION OF ONSAGER 
COEFFICIENTS 

After the preliminaries in the previous sections, we can 
now calculate the Onsager coefficients for our finite-time 
Carnot cycle in the linear-response regime AT — > as 
follows. The first step is an appropriate choice of the 
thermodynamic forces and fluxes for this system. A typ- 
ical way to choose these thermodynamic forces and fluxes 
is to introduce the rate of the total entropy production 
a during one cycle: 



Qh Qc 

~T h ~ T c 



(11) 



This is just the entropy increase of the two reservoirs 
during one cycle divided by the cycle period 2(1% — l\)/u 
because the entropy of the working substance does not 
change after one cycle. Using the relation Q c = W — Qh 
and considering the linear-response regime AT — * 0, it 
can be rewritten as 



u{-W) 
2(l 3 -h)T 



AT 

Ji2 



(12) 



where T = (T h +T c ) /2 and we have neglected AT 3 Q h and 
uWAT terms, the reason of which we will clarify later. 
According to the linear-response theory, & can be ex- 
pressed as the sum of the product of the thermodynamic 
force and its conjugate thermodynamic flux (TH. fl6j| : 



a = J\X\ + J2X2, 
where we define the thermodynamic forces as 
—W 



2{h-h)r 



AT 

J>2 



(13) 



(14) 



Ji = u, J 2 



(15) 



Moreover the linear-response theory assumes the Onsager 
relations between the fluxes and forces [3, [T3| : 



Ji 



u = L\ 

Qh = -^21 



-W 
~W 



AT 

Luj^, (16) 



AT 

2(l 3 -h)T + 122 T*> (17) 



where Ly 's are the Onsager coefficients and the non- 
diagonal elements should satisfy the symmetry relation 
L12 = L2i- From these relations between the thermody- 
namic fluxes and forces, we understand that & — J\X\ + 
J2X2 is the quantity of the second-order of the ther- 
modynamic forces, which explain why we neglected the 
higher order terms like AT 3 Q h and uW AT in Eq. 1(12]). 
Moreover, although we considered the contributions of 
the additional heat transfers QjJ dd and Qc dd to the to- 
tal heat and the work in Sec. Ill, we can easily show 
that their effects do not contribute to a in the limit of 
AT — > 0. Therefore we can indeed neglect them in the 
linear-response regime and use Qh = Qh an d W — W st 
in the calculations of the Onsager coefficients below. 

We are now in a position to calculate the Onsager coef- 
ficients explicitly. First we determine L\\ and L21 as fol- 
lows. To calculate Tn, we consider the relation between 
u and X\ in the case of AT = 0. Expanding Eq. ([9]) by 
AT and putting AT = 0, we can obtain the relation 



W = —AmNA\ 



1 L \ V 2 



(18) 

Comparing Eq. (fTS)) with Eq. (Tl6|) . L u is determined as 
(is - h)T 112 



lrnN(l + ±)J 2 -^ln% 



(19) 



Likewise Qh with AT = can be evaluated up to the 
linear order in W using Eq. © and (fTS)) as 



Qh = 



fc R T 3 / 2 



-W 



4m(i + ^) 



27rfe E 



2{h-h)T 



(20) 



Comparing Eq. ([20| with Eq. (|17ll . L21 is determined as 



fc R T 3 / 2 



4 m (I + l ) /i^fe 



(21) 



Next we determine T12 and T22 as follows. To calcu- 
late Ti2, we consider the relation between u and X2 in 
the case of W = 0. In this case, regardless of a finite 
temperature difference, useful work cannot be obtained 
because the engine runs so fast that it cannot output 
positive work. Therefore, this case is called the work- 
consuming state. The speed of the piston at the work- 
consuming state can be obtained as a solution of W = 
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in Eq. 



(|9|). Considering in the linear-response regime 
V% ~ AT/(2y/T), it becomes 



fc R T 3 / 2 



AT 



M£ + &) 



27rA: E 



T 2 



(22) 



From Eq. (|16j) and Eq. ([22]), we can obtain 

fc R T 3 / 2 



^12 — 



M£ + &) 



27rfc E 



(23) 



From Eqs. ([21]) and ([23]) . we can confirm the symmetry 
relation L±2 = T21 as expected. To determine the last 
coefficient T22, we consider the heat flow Qh at the work- 
consuming state using Eq. (JSJ) and 



Qh = 



TVfcllnf 



8m (i 



_j.5/2 



AT 

y2 ' 



(24) 



From Eq. (fTF]) and Eq. (|24|) . the last coefheient L22 turns 
out to be 



£22 — 



Nkl In ^T 5 / 2 



(25) 
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FIG. 3: T-dependence of the Onsager coefficients. The dashed 
line, the solid line and the dotted line indicate the theoretical 
Onsager coefficient L 22 (Eq. (J25J)), L 12 (Eq. J23J) (= L 2 i 
(Eq. (HlJ)) and in (Eq. (]T9)l '). respectively. The parameters 
used in the MD simulations arc AT = 4 x 10" 3 and d = 0.01. 
The other parameters are the same as in Fig. [2] The MD data 
were obtained by averaging 100-800 cycles after transient 5 
cycles in the simulations. 



Note that the positivity of the rate of the total entropy 
production, & = J\X\ + J2X2 > 0, should restrict the 
values of the Onsager coefficients to L\\ > 0, L22 > 
and L11L22 — L12L21 > 0. We can confirm that of 
our finite-time Carnot cycle surely satisfy these relations. 
These analytic expressions of the Onsager coefficients Lij 
are the main result of this paper. 

To confirm the validity of the above analytic calcu- 
lations of the Onsager coefficients, especially T depen- 
dence of them, we performed the event-driven molecular 
dynamics (MD) computer simulations 0, [H, [13] of our 
two-dimensional finite-time Carnot cycle, following the 
procedure described in Sec. II. 

To calculate L12 and L22 at given T, we fix the tem- 
perature difference AT to a sufficiently small value and 
find the piston speed where the work becomes 0. Then 
we can determine T12 and T22 numerically as 



T12 
T22 



Jx _ uT 2 
~X~2 ~ AT' 

h_ = QhT 2 

x 2 



AT 



(26) 
(27) 



from Eq. (flB) and Eq. (|17|) . Next, to calculate L\\ and 
Li2\ at given T, we set AT = 0. Fixing u to a sufficiently 
small value (2.5 x 1CT 5 - 3 x 1CT 4 ), we can determine 
L11 and L21 numerically as 



^11 



L, 



2(^3 -h)Tu 
W 

2(l 3 - h)TQ h 
W 



(28) 
(29) 



from using Eq. (|16p and Eq. (jTTJ) . Fig. [3] shows T depen- 
dence of these Onsager coefficients determined by the MD 
simulations as well as the analytic results. We can see 
fairly good agreement between the MD data and the an- 
alytic lines in the range we studied. These data clearly 
support the validity of our analytic expressions of the 
Onsager coefficients which has been obtained under some 
theoretical assumptions. 



THE EFFICIENCY AT THE MAXIMAL 
POWER 



To see how the Onsager coefficients derived in Sec. IV 
indeed govern the behavior of the finite-time Carnot cy- 
cle, we briefly introduce the general framework of the 
heat engine obeying the Onsager relations Q. 

The basic setup is as follows. We consider a general 
steady or cyclic process in which the work is extracted 
from the heat flow between the small temperature differ- 
ence, (see Fig. |4|) The work W done against the external 
force F is W — —Fx where x is the thermodynamically 
conjugate variable of F. We define a thermodynamic 
force as X\ = F/T and the corresponding thermody- 
namic flux as Ji = x. We also choose X2 — 1/T C — 1/Th 
as another thermodynamic force and J2 = Qh as the 
corresponding thermodynamic flux. If we consider the 
linear-response regime AT — > 0, X 2 can be written as 
X2 — AT/T 2 . Moreover the linear-response theory as- 
sumes the Onsager relations between the thermodynamic 
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Qi 



T h =T+^f 



T r = T-^ 



4 



W 



Qc 



FIG. 4: Schematic illustration of the heat engine governed by 
the Onsager relations Eq. (|30p and Eq. (|31[l . 



forces and fluxes [HI, : 



Ji — LuXi + L12X2, 
Ji — L21X1 + L22X2, 



(30) 
(31) 



where the non-diagonal elements of the Onsager coeffi- 
cients should satisfy L 12 = L 2 i- Then, the power P = W 
and the efficiency 77 = P/Qh of the engine can be ex- 
pressed as 



P 
V 



-JiXtT, 
-J1X1T 
J2 ' 



(32) 
(33) 



Now the efficiency at the maximal power ?7 max can be 
given as follows: We first maximize the power at X\ = 
— Li2^2/(2I/n) which is determined as the solution of 
dP/dXi = 0, then, ?7 max becomes 



AT 



?7max 



q- 



where 



2T 2 



-^12 

^/LllL22 



(34) 



(35) 



is called the coupling strength parameter. Note that it 
takes — 1 < q < +1 due to the positivity of a. When q 
satisfies the tight-coupling condition |g| = 1, ?7 max takes 
the maximal value AT/(2T), which is equal to the CA 
efficiency up to the lowest order in AT. 

In fiefs. 0, [f|, we studied ry max of this finite-time 
Carnot cycle extensively by performing the MD computer 
simulations as a numerical experiment to verify the valid- 
ity of the CA efficiency. We found there that ?7 max agrees 
with the CA efficiency in the limit of AT — * 0. Our 
molecular kinetic theory also confirmed this property Q : 
Using Qi — Qf + Qf dd in Sec. Ill, we can calculate the 
speed u = u max at which the power P maximizes as 



k B (T h -T c ) 



4/ " /l ' i + \ 




8Lk B 
3SAV% 

§or*4 (36) 



and can confirm that the efficiency 77 at u 



shows 



l(Hmax) = 



W(u max ) 



Q\i ( M max) Qfo (^max) 



AT 
~2T 



(AT-> 0). 



(37) 



Now we can clarify the underlying physics of this be- 
havior of ?7 max . In the limit of AT — * 0, our finite-time 
Carnot cycle can be described by the Onsager relations 
as shown in Sec IV. Then we can confirm that q = 1 in 
our finite-time Carnot cycle from Eqs. (fTSj) . (|2"3")l . 
and (|35p. This condition gives a proof that our finite- 
time Carnot cycle shows the CA efficiency in the limit of 
AT — * as suggested in [3, [E[ from the viewpoint of the 
linear-response theory. 



VI. SUMMARY 



In this paper, we have studied a finite-time Carnot cy- 
cle of a two-dimensional weakly interacting nearly ideal 
gas working in the linear-response regime and have ex- 
plicitly calculated the Onsager coefficients of this sys- 
tem for the first time. Molecular dynamics computer 
simulations of this cycle have supported the theoret- 
ical calculations in spite of some assumptions in the 
analysis. We have revealed that the Onsager coeffi- 
cients of this system satisfy the tight-coupling condition 
q = L12/ v ' L11L22 — 1 and therefore can understand why 
our finite-time Carnot cycle attains the Curzon-Ahlborn 
efficiency in the linear-response regime AT — * as sug- 
gested in Refs. [3, It would be an interesting prob- 
lem to construct and study the collective behavior of 
the Carnot cycles coupled with each other by using the 
property of the Onsager coefficients derived in this pa- 
per @, EE m, [10, EH . 
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